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Abstract 

We re-derive the renormalization group equation for the effective couphng of the 
dimension five operator which corresponds to a Majorana mass matrix for the Stan- 
dard Model neutrinos. We find a result which diff^ers somewhat from earlier calcu- 
lations, leading to modifications in the evolution of leptonic mixing angles and CP 
phases. We also present a general method for calculating /^-functions from counter- 
terms in MS-like renormalization schemes, which works for tensorial quantities. 
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1 Introduction 



The Standard Model (SM) is most likely an effective theory up to some scale 
A, above which new physics has to be taken into account. The discovery of neu- 
trino masses requires an extension of the SM, which may involve right-handed 
neutrinos or other new fields. Introducing right-handed neutrinos allows Dirac 
masses via Yukawa couplings analogous to the quark sector. In general, 
lepton number need not be conserved, so that Majorana masses are possi- 
ble. For left-handed neutrinos this can, for example, be achieved with Higgs 
triplets. Right-handed neutrinos can have explicit Majorana masses Mr of 
order A, since they are gauge singlets and since there are no protective sym- 
metries. This leads to a picture with zero or tiny left-handed Majorana masses 
Ml, with mj) similar to the charged lepton masses, and with a huge Mr. Di- 
agonalization of the neutrino mass matrix results in Majorana fermions and 
eigenvalues ~ Mr, and Ml — m^/AfR. For Ml = the neutrino masses are 
thus given by the see-saw relation m^/MR [1], which provides a convincing 
explanation for the smallness of neutrino masses. 

Another, less model dependent approach is to study the effective field the- 
ory with higher dimensional operators of SM fields. If lepton number is not 
conserved, some of these generate Majorana neutrino masses. The lowest di- 
mensional operator of this kind has dimension 5 and couples two lepton and 
two Higgs doublets. It appears e.g. in the see-saw mechanism by integrating 
out the heavy right-handed neutrinos. 

As quarks have only small mixings, it is somewhat surprising that neutrinos 
most likely have two large mixing angles [2-4]. It is interesting to investigate 
mechanisms which can produce such large or maximal mixings. These mecha- 
nisms operate, however, typically at the embedding scale A. For a comparison 
of experimental results with high energy predictions from unified theories, it is 
thus essential to evolve the predictions to low energies with the relevant renor- 
malization group equations (RGE's). This evolution is related to the running 
of the leading dimension 5 operator. Therefore, we calculate in this letter the 
RGE that governs this running above the electroweak scale at one-loop order 
in the SM. 



2 Lagrangian and Counterterms 

Let ii, f e {1,2,3}, be the SU(2)L-doublets of SM leptons, the SU(2)l- 
singlet (right-handed) charged leptons, and 4> the Higgs doublet. The dimen- 
sion 5 operator that gives Majorana masses to the SM neutrinos is given by 
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(1) 



where k is symmetric under interchange of the generation indices / and ^f, £ is 
the totally antisymmetric tensor in 2 dimensions, and '■— (^l)*^ is the charge 
conjugate of the lepton doublet. a,b,c,d e {1,2} are SU(2) indices. They 
will only be written explicitly in terms with a non-trivial SU(2) structure. 
Summation over repeated indices is implied throughout this letter. 

gives rise to the vertex shown in Fig. 1, and an analogous one for the 
Hermitian conjugate term. 



Fig. 1. Vertex from the effective dimension 5 operator and the corresponding Feyn- 
man rule. The gray arrow indicates the fermion flow as deflned in [5]. 

The complete Lagrangian consists of the SM Lagrangian =$fsM and proper 
counterterms 



In the following, we omit most of those parts that yield only flavour diagonal 
contributions to the /3-function and therefore do not contribute to the running 
of mixing angles, in particular terms involving quarks and gauge bosons. The 
remaining ones are 




(2) 




-(i;We^V^{ + h.c. ; 



(3a) 
(3b) 
(3c) 



^kin(€0=4'(^7^a^)(5^^J,X, 

^Yukawa = -((5Fe)g/ eS^^f^ + h.C. ; 



(4a) 
(4b) 
(4c) 



+ h.C. . 



(5) 
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5Zi {i G {^Lj 0}) determine the wavefunction renormalization constants Zi — 
1 + 5Zi, defined in the usual way. Note that Z(,^ is a matrix in flavour space. 
5k satisfies the relation 

«B = Z--- {Ziy'' [k + Sk] z;}z-y- , (6) 

where the factor /x^ is due to dimensional regularization, with ^ denoting the 
renormalization scale and e := 4 — The subscript B denotes a bare quantity. 
Note that the usual ansatz kb ~ Z^^^k is not possible in this case, as it would 
obviously spoil the symmetry of or k with respect to interchange of the 
flavour indices. 



3 Calculation of the Counterterms 



In the MS scheme, the quantity 5k can be computed at one-loop order from 
the requirement that the sum of diagrams in Fig. 2 be ultraviolet finite. 



" +f m K + I — 





+ diagrams from gauge interactions 



Fig. 2. Diagrams relevant for the renormalization of the vertex from the effective 
dimension 5 neutrino mass operator. The last diagram represents the counterterm. 



Using FeynCalc [6] we obtain 



(7) 
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where denotes the contribution from gauge interactions. The usual calcu- 
lation of the wavefunction renormalization constants yields 



(8) 
(9) 



Again, and C^^ represent terms from quarks and gauge interactions, which 
are diagonal in flavour space. 



4 Calculating RGE's from Counterterms with Tensorial Structure 



The calculation of the /5-function involves some subtle points, which are related 
to the matrix structure of the counterterm Lagrangian. Before presenting our 
result in Sec. 5, we provide now some details of the calculation, which should 
be of general interest and which are essential for verifying our result. In partic- 
ular, we generalize the usual formalism for calculating /9-functions to include 
tensorial quantities as well as non-multiplicative renormalization. 



We are interested in the /9-function for a quantity Q, /5q := A*-^- In general, 
the bare and the renormalized quantity are related by 



Qb = z^i ■■■z;^[Q + 5Q]^i^'^^ z;^^i . . . z;- 




n^;; [^3+5^]/^^'=^ [UKi] ' (10) 

Kiel ) 



where / = {!,..., M}, J = {M -|- 1, . . . , A^} and Dq is related to the mass 
dimension of Q. 5Q and the wavefunction renormalization constants depend 
on Q and some additional variables {Va}, 



5Q = 5Q{Q,{Va}), (11a) 
Z^^ = Z^^{Q,{Va}) {l<i<N). (lib) 

Note that Q = Q{[x) and Va = Va{ij) are functions of the renormalization scale 
/X, but 5Q and Z^^ do not depend explicitly on fi in an MS-like renormalization 
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scheme. Taking the derivative of equation (10) yields 



I ,_,-I^Q^^±Q^ 

an 



\iei 



/d6Q 
\ dQ 



+ 



+E 



'd6Q 



dVA 



(3y\+eDQ{Q + 5Q) 



Vie/ / I ieJ \i'<i 



X 



X 



'dz: 



dQ 



jeJ \j'<j 

dz: 



E 



dV. 



\j">j 



iel \i'<i 



'dz: 



dQ 



+ E 



'dZ', 



dVA 



Va 



V'>i J ) ViGJ / 



X 



(12) 



Here we have introduced the notation 



'dF 
dx 



dF 

dx 



y 



for scalars x, y 

for vectors x = {xm),y = {Vm) 



^ dF 

dF 

E 1 y^nn for matrices x = {xmn),y = {ymn) 

m,n '^^mn 

etc. . 



(13) 



We will solve equation (12) and the corresponding expression for Va by ex- 
panding all quantities in powers of e. In the MS-scheme the quantities 5Q and 
Zj,. can be expanded as 



5q = j: 



k>l 



6Z. 



(14a) 
(14b) 



k>l 



with higher powers of ^ corresponding to higher powers in perturbation theory. 
On the other hand, /3-f unctions are finite as e ^ 0. We can therefore make 
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the ansatz 



(15a) 
(15b) 



where n is an arbitrary integer. Note that in this case the power of e is not 
related to the order of perturbation theory. Prom (14) and (15) we find that 



dQ 



-idZ4 



dg 



d6Z^ 
dQ 



(16) 



where the lowest possible power of ^ appearing on the right side of (16) is 1. An 
analogous relation holds for Q ^ Va- Our analysis of equation (12), starting 
with the inspection of the e" term, then shows that f3Q^ vanishes. The analog 

of equation (12) for /5y^ implies that f3y^ vanishes as well. Repeating this 
argument for successively smaller positive powers of e implies that 



= -sDqQ , 



yke{2,...,n}, 



PVa 



-eDv^VA 



(17a) 
(17b) 
(17c) 



Note that these terms do not contribute to the /^-function in 4 dimensions, i.e. 
for e — > 0, but they are necessary to read off (3q^ from equation (12), leading 
to the result 



(0) 



Q 



/d6Q^ 
\ dQ 



+ E^i 

is/ 



D. 



Q^ 

Do 



A 

'dZ, 



dQ 



\ dVA 
Q) + Y.Dv, 



Va)-Dq5Q,, 



/dZ^^,i 




Va^ 


\ dVA 



'dZ„ 



Q 



dQ 



q)+j:d 



dZ, 



Va 



dVA 



Va 



■Q 



(18) 



Note that for complex quantities Q and Va we have to treat the complex 
conjugates Q* and VX as additional independent variables. 



5 Renormalization Group Equation 



The RGE for the effective coupling k, is 

.^ = &. (19) 
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Using equations (18) and (7) - (9), wc obtain for the contributions from vertex 
and wavefunction renormalization (omitting terms from C^, and C^l)- 



+ \k 



igtt^/?!^) = -2 [k{y}y,) + (y;ty;)' 



(20a) 
(20b) 



Adding tlie terms involving quarks and gauge bosons [7,8], we obtain tiie final 
result 



+ 



+\k - zgln + 2 iv (syj-y, + W^Yd + Y}Y,) h , 



(21) 



where g2 is the SU(2) gauge coupling constant and where Y^, Y^ are the Yukawa 
matrices for the up and the down quarks. Thus, compared to earlier results 
[7], we find a coefficient — | instead of — | in front of the non-diagonal term 
K(Y^Ye)-\-{Y^YeY K. Note that the difference in the Are-term is due to a different 
convention for the Higgs self-interaction used in this work. 

We have checked our results by calculating the essential parts of the same (3- 
functions from the finite parts of the relevant diagrams in the framework of an 
underlying renormalizable theory. This calculation as well as the application 
to the MSSM and the two Higgs SM will be presented in a future paper [9] . 



6 Discussion and Conclusions 



We have calculated in the SM the /3-function for the effective coupling k of 
the dimension 5 operator which corresponds to a Majorana mass matrix for 
neutrinos. We have explicitly presented our calculations for the non-diagonal 
part of the /3-function, where our result disagrees with the previous one in [7] 
by a factor of 3. This part is responsible for the evolution of neutrino mix- 
ing angles and CP phases. Therefore, our result modifies the renormalization 
group running of these quantities between predictions of models at high en- 
ergies and experimental data at low energies. Consequently, our work affects 
the SM results of previous studies based on the existing RGE's, e.g. [10-15]. 
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